Deformation of Quantum Dots in the Coulomb Blockade Regime 
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We extend the theory of Coulomb blockade oscillations to quantum dots which are deformed by 
the confining potential. We show that shape deformations can generate sequences of conductance 
resonances which carry the same internal wavefunction. This fact may cause strong correlations of 
neighboring conductance peaks. We demonstrate the relevance of our results for the interpretation 
of recent experiments on semiconductor quantum dots. 
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Recently two experimental groups have reported the 
first measurements of the distribution of Coulomb block- 
ade peak heights in semiconductor quantum dots [0,0. 
These measurements were designed to test various pre- 
dictions of the statistical theory of the Coulomb block- 
ade [pi. This theory was formulated by Jalabert, Stone 
and Alhassid for the low temperature regime kT < A, 
where A is the mean resonance spacing. In this regime, 
transport through the quantum dot occurs by resonant 
tunneling through a single quasibound state. On the ba- 
sis of a random matrix assumption for the quasibound 
eigenf unctions of the dot, the statistical distribution of 
the peak heights of the conductance resonances was de- 
rived yj both for systems with time-reversal symmetry 
(vanishing magnetic field) and for systems with broken 
time-reversal symmetry (with sufficiently large magnetic 
field). 

In the experiments [|l|,|[, good agreement was found 
with the predicted form of the distribution functions. At 
the same time, unexpected correlations were observed g] 
between neighboring peak heights. The correlations com- 
prised up to four adjacent peaks. Even stronger cor- 
relations have recently been found in a set of experi- 
ments [|| addressing the phase of the transmission am- 
plitude through a quantum dot. In these experiments, 
more than 10 consecutive conductance peaks displayed 
the same phase, and nearly the same peak height. This 
phenomenon has not been explained so far. 

In this paper we propose a novel mechanism which can 
account for the above-mentioned correlations. It applies 
whenever the same set of capacitively coupled gate elec- 
trodes is used to define the depth and shape of the dot 
confining potential. This is the case in the experiments 
of Refs. [|l|,||,0 but not for others using a back gate. We 
study the effect of a deformation of the potential on the 
eigenvalues and eigenfunctions of the dot. 

Several authors B have shown that a parabolic poten- 
tial is a useful approximation to the confining potential 
of a quantum dot. To simplify the argument, we start 
with this case. We identify sequences of conductance 



resonances which, as a result of level crossings, carry the 
same internal wave function. Similar sequences are found 
if the parabolic potential is replaced by a more realistic 
confinement potential with level repulsion and chaotic 
classical electron motion. We first disregard the charging 
energy of the dot but show later that our conclusions re- 
main essentially unchanged when this energy is included. 
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FIG. 1. Energy levels of the two-dimensional harmonic os- 
cillator with frequencies cux ~ 1 and u>y — 1 — P as functions 
of/?. 



We choose the x axis through the center and in the 
direction of transmission through the dot. To account 
for anisotropics imposed by the external gate potential 
(gp) , we choose different values for the frequencies lUx , ^y 
of the parabolic confining potential. For simplicity, we 
keep the frequency u>x = w^ fixed and change the fre- 
quency u!y = uj^{l — (3), with Q < (3 < 1. Increasing /3 
widens the potential in the j/-direction and thus deforms 
it. Fig. 1 shows part of the single-particle spectrum ver- 
sus f3. The shell structure at /? = is clearly visible. 
The energy levels form a network of intersecting straight 
lines. Some lines run nearly parallel to the /3-axis while 
others have large negative slopes. The explanation is 
simple. Each single-particle state is characterized by two 
quantum numbers (rix, ny). States with small (large) val- 
ues of Uy have small (large) negative slopes. When two 
lines cross, the quantum numbers are carried along the 



straight lines. We show below that the pattern in Fig. 1 is 
generic and applies also to the single-particle spectrum of 
deformed non-integrable potentials, with one difference. 
The points of intersection of the straight lines disappear 
and are replaced by avoided (Landau-Zener) crossings. 

To present the essence of our argument, we replace 
Fig. 1 by an idealized picture (Fig. 2): A set of equally 
spaced straight lines (A-levels) runs parallel to the /3-axis, 
a second similar set (_B-levels) has large negative slope. 
To model the generic case, we have replaced actual cross- 
ings by avoided crossings. Single-particle wave functions 
retain their identity across avoided crossings in nearly 
the same way as if the lines were to intersect. Therefore, 
the wave function on any A-level is nearly independent 
of /3, and so is the Fermi energy Ep marked by a dashed 
line. At a fixed value of /3, all single-particle states up 
to Ep are filled. A change of the gp (and, therefore, 
of (3) is slow on the scale of characteristic times of the 
quantum dot. Therefore, electrons in occupied orbitals 
follow the deformation (3 adiabatically. Let 5(3 be the dis- 
tance between i?-lines, and g'psfi be the change of the gp 
needed to change j3 into (3 -f 5(3. What happens as we in- 
crease the gp by gP5/3? The last occupied A- level moves 
adiabatically one "floor" down to its nearest analogue, 
thereby changing its wave function. The A-level right 
below Ep becomes empty. If gp5/3 =gPe, the change of 
the gp needed to pull another electron into the dot, then 
the next conductance resonance will carry the very same 
single-particle state as its predecessor. (The Figure sim- 
plistically suggests that the next resonance occurs when- 
ever a solid line crosses Ep] in reality the criteria given 
in gl apply). This process can repeat itself for a second, 
third, etc. time. Thus, subsequent conductance peaks 
might not be independent, but would be manifestations 
of one and the same single-particle state. 
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FIG. 2. Schematic illustration of the mechanism described 
in the text. The arrows indicate the levels which become avail- 
able for occupation when the deformation is increased. 

This picture raises several questions, and we devote 
the remainder of this paper to the answers, (i) How real- 
istic is Fig. 2? (ii) Does a significant deformation indeed 
occur, and is it reasonable to assume g'psp '^ gPe? (iii) 
How does the charging energy affect our mechanism? 



(i) To construct a more realistic case without the reg- 
ularity of Figs. 1 and 2, we consider the two-dimensional 
single-particle Hamiltonian 

H = ^ + ^ix' + il- P)V) )^ncoL', (1) 
2m 2 

where L = iiy-^—x-^) is the dimensionless z-component 
of the angular momentum operator. In nuclear physics, 
the three-dimensional analogue of H is known as the Nils- 
son model and has been quite successful in explaining the 
spectra of deformed nuclei 0. For /? > and X ^ 0, H 
is not integrable ||§] and displays level repulsion. We may 
view iJ as a mean-field approximation to the total Hamil- 
tonian. The latter contains both, the mutual interaction 
of the electrons and the real confining potential. Any 
such approximation should lead to more or less chaotic 
single-particle motion M. 
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FIG. 3. Energy levels as a function of /3 of the Hamiltonian 
(0) for A = 0.004. 

Fig. 3 shows part of the spectrum of H versus /3. For 
/3 = the eigenvalues m oi L are good quantum numbers. 
The levels with m ^ are pairwise degenerate. They 
form ascending sequences with ±m, ±(m — 2), ... , ±1 or 
0. For /? > 0, each eigenfunction contains a superpo- 
sition of m-values. Depending on the strength A, the 
mixing may extend over m-values from different shells. 
The overall pattern is quite similar to Fig. 1. The resolu- 
tion in Fig. 3 is not fine enough to show that all crossings 
are avoided crossings. But is it realistic to assume that 
wave functions on (nearly) horizontal levels retain their 
identity when (S is changed? This, after all, is the cen- 
tral point used in our argument. We have inspected the 
eigenvectors for the case A = 0.005 in the basis defined 
by the eigenvectors of the undeformed Hamiltonian. The 
levels with the smallest slope have a dominant compo- 
nent. It is the state with quantum numbers ny — and 
Ux = A'shcn- Here, iVshcii denotes the shell number. The 
modulus of the amplitude of this component is 0.85 or 
larger and remains virtually unchanged for (S < 0.2. Ad- 
mixtures of non-dominant components with amplitude 
moduli greater than 0.1 are limited to a small number of 
states (ten or less). In contradistinction, the eigenvectors 
of the steepest levels are genuine mixtures without dom- 



inant components. For /? > 0.1, the number of essential 
components (amplitude moduli larger than 0.1) may be 
close to one hundred. There is a gradual transition, of 
course, from the levels with the smallest to those with 
the largest slope. 

The wave functions of the dominant components cor- 
respond classically to orbits which oscillate in the x- 
direction, with little or no motion in the y-direction. The 
relative stability of these orbits is reflected in the small 
amount of mixing of their quantum analogues. We ex- 
pect the bulk of the current to be carried by such orbits. 
The experimental set-up appears to enforce occupation 
of the associated characteristic levels as the conducting 
agent; otherwise a current could barely flow. 
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FIG. 4. Sections of the confining potential in (a) the 
a;-direction and (b) the y-direction in a plane parallel to the 
surface at a distance of 0.01 ^m. The rods on the surface 
are arranged as shown in the inset. The distance from top 
to bottom rod is 2 /xm, the drawing is in scale. The three 
curves in (a) and (b) refer to a plunger potential (middle left 
rod) of —!/(), —0.9Vo and — 0.8Vb, respectively, where — Vo is 
the potential at the other rods. 

(ii) Is it reasonable to discuss deformation, and how 
does variation of the plunger gate voltage translate into 
a change of (31 A quantitative answer would require pre- 
cise knowledge of a particular experimental set-up, and 
we must confine ourselves to semiquantitative considera- 
tions. We consider a set of gates arranged as indicated in 
the inset of Fig. 4b. The current flows in the vertical di- 
rection. The middle left conducting rod serves as plunger 



gate. To study the potential well V{x, y) formed beneath 
the gates, we have taken a set of line charges. This avoids 
solving the full boundary-value problem. The charges on 
the top, middle right and bottom rods were all chosen 
equal. Upon changing the charge on the plunger, they 
were re-adjusted to ensure that the potential values at 
the height of the entrance and exit barriers and at the 
tip of the middle right rod remain unchanged. An in- 
crease of the plunger voltage (charge) by 20% results in 
a decrease of 5% and 0.7% at the left hand and the right 
hand rods, respectively. Sections of the potential well 
confining the electrons taken in a plane below the gates 
are shown in Fig. 4, together with the changes due to an 
increase of the plunger voltage. In the x-direction, the 
bottom of the potential is lowered. In the y-direction, 
the same lowering is accompanied by a widening toward 
the left of the potential. Figure 4 indicates clearly the 
importance of deformation. 

To describe these changes quantitatively, we approxi- 
mate the potential by a paraboloid. The lowering of the 
bottom results in an effective increase of uj^. In the y- 
direction, there is an effective decrease of ujy. To first 
order in a, we quantify these changes as tOx — > (1 -I- q;)^^; 
and LOy ^ {1 — ct)^y Rescaling the energy by a factor 
(1 -I- a) and writing /3 = 2a, we obtain Wj, — > (1 — l3)ujy 
while LUx is fixed. This is the situation discussed above. 
With — Vb < the original value of the voltage on all the 
gates, an increase of the plunger potential from —Vq to 
— 0.75Vb (a realistic range of values) causes, in the exam- 
ple in Fig. 4, a deformation f3 « 0.1. This is a sizeable 
value and corresponds to the scales shown in Figs. 1 and 
3. 

To relate gps/s and gpe, we estimate the change of 
depth of the harmonic oscillator potential needed to pull 
another electron into the dot. As shown in Fig. 4a, we 
require that the values of y(±a;o,0) remain unchanged 
under changes of the plunger voltage. Here, ±xq denote 
the X- values of the positions of the barriers which define 
entrance and exit of the dot. We take y(±xo,0) = Ep- 
(Since there are barriers, the value of Ep is actually 
slightly lower, but this is immaterial for what follows). 
In the case of Coulomb blockade without deformation [^ , 
addition of one electron requires the potential to be low- 
ered by As+e'^/C « e^/C. Here, A^ is the average single- 
particle level spacing and e^/C the charging energy. For 
N electrons, this yields Ne^/C [|lO[. Measured from 
the points ±a;o, the parabolic potential in x-direction 
must therefore have the depth (1 + a)'^{ujxXo)'^/2 « 
Ep + Ne'^/C. Inclusion of deformation modifies the value 
of C but does not affect this argument qualitatively. We 
combine this result with the corresponding expression for 
the case of iV -I- 1 electrons at deformation a + 1/25(3 to 
obtain 5(3 ~ 2{e'^ /C)/{ujxXi^Y ^ l/N. In the last step we 
have neglected Ep in comparison with Ne^ /C . We com- 
pare this result with the average distance in (3 of nearest 
avoided crossings. We consider the number L of cross- 



ings on a horizontal level at A'shcii in Fig. 1. We find 
that L is governed by [/?/(! - I3)]N'^^^J2. For the dis- 
tance of crossings this yields 5(3 ^ 1/^shcii ^ 1/^- This 
result is generic. The estimates show that roughly one 
level crossing occurs for each additional electron pulled 
into the dot. This result lends strong credibility to our 
model. 

(iii) Figs. 1 to 3 represent the single-particle energies 
without taking account of the charging energy Uq — 
e^/C. To redress this omission, we use the framework 
of the self-consistent treatment of the Coulomb interac- 
tion within the quantum dot as applied in Ref. [0|. In 
the ordinary mechanism of avoided crossings, two levels 
tjiy) with J = 1, 2 are functions of an external parame- 
ter V . Let the two levels be coupled by a (real) matrix 
element W . The eigenvalues are 



E^='-l±^± 



ei -£2 



W^. 



(2) 



If the two levels ^j{V) intersect at some value V = Vq, 
then the eigenvalues E± show an avoided crossing &tV = 
Vq. We now include the charging energy by considering 
two "effective" single-particle levels E± defined as solu- 
tions of the mean-field equations lO] E± = E±-\'Uo{n^) . 
Here the {n±) are determined self-consistently. 



(n±) = - 

TT 



^ (2{Ef-E±)\ tt 
arctan 



(3) 



where T± are the intrinsic widths of the two levels. We 
solve these equations sclf-consistently for E± and {n±) 
and pick the solution that minimizes the mean-field en- 
ergy Emf ~ {n+)E+ + {n^)E^ - q{n+){n^). We focus 
on a situation where the two effective levels are within an 
interval C/q around the Fermi energy, with one level below 
Ep and the other above it. As a result of the charging en- 
ergy, the two levels participating in the avoided crossing 
are simply shifted apart by roughly C/g- The wave func- 
tion at y <C Vb of the lower (higher) state is the same as 
the wave function at T^ ^ Vb of the higher (lower) state. 
The mixing of the two wave functions near the avoided 
crossing is independent of the charging energy Uq. As 
a result, the wave function of the state that is occupied 
for V^ <^ Vq changes smoothly into the wave function of 
the state which is empty, and available for the next elec- 
tron, for V ^ Vq. This is in keeping with the mechanism 
proposed above. 

Clearly, strong correlations must be expected for the 
conductance peak heights and for the transmission phases 



within a sequence of resonances caused by one and the 
same quantum state. Long sequences of this type are 
more likely to occur in quantum dots with regular shape. 
This might explain why long sequences of similar reso- 
nances where observed for the symmetric (regular) dot 
of Ref. B and perhaps weaker correlations in the more 
irregular dot of Ref. 0. Our results suggest that fur- 
ther correlations not tested so far in Ref. B might exist. 
This possibility as well as a study of the effects of an ap- 
plied magnetic field will be investigated in a forthcoming 
paper. 
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